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ABSTRACT 


A simple analytic model is formulated which should prove useful for 
evaluating a satellite surveillance system. This model describes the 
fraction of orbits in which a long lived satellite will "see" a selected 
position on the earth's surface. The fraction of orbits, denoted by P, 
is a function of the latitude of the ground position in question, the 
orbital inclination of the satellite, and the sweep width of the satel- 
lite's sensor. 

P is used to find estimates for the mean time to next sighting of a 
selected ground position given that position has suddenly become of in- 
terest, to estimate the number and locations of ground readout stations, 
and to estimate the number of times a satellite would sight a transitting 
object on the earth's surface. In addition, the application of P to the 
problems of optimum orbit, evasive countermeasures, and satellite weapon 


delivery is discussed. 
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low Intmoduction. 
Rapid developments in space technology have resulted in.an increased 


use of unmanned earth-orbitting satellites for both civilian and military 





purposes. These satellites provide information on outer space, the earth's 
atmosphere, and the earth's surface. Satellites for earth surveillence 
offer a distinct advantage over earth-bound systems because of their abil-~ 
ity to scan large segments of the earth in a relatively short period of 
time, and because they are unencumbered by national air-space restrictions. 

The Tiros and Nimbus series weather satellites represent well known 
examples of the use of satellites by the United States for earth surveil- 
lence. Tiros III, launched in 1961 was equipped with photographic sensors 
and provided information on storms and storm paths. Nimbus I, launched in 
1965 served a similar purpose but made use of television cameras as sen- 
sors. 

Although weather received the first attention, it is not the only 
phenomenon worthy of satellite surveillence. Nuclear test explosions are 
of sufficient interest that two Vela Hotel satellites were launched by 
the United States in 1965. Additional satellites apparently have been 
launched for certain classified military purposes. The Samos series 
satellites evidently have, for the last five years, been photographing 
military installations in Russia and Red China. [6] Seven military recon- 
naissance satellites (possibly of the Samos series) were reported to have 
been launched during 1965. [4] These satellites had polar orbits and 
mission lives of about five days. The observations made by these satel- 
lites are unknown, but could include ship and troop movements or test 
firings of ICBM's, as well as military installations. 


Given a surveillence task, alternative satellite systems might be 


designed to perform this task. We may want to decide which system or com- 
bination of systems is most cost-effective. Or we may want to compare 

the alternatives with more conventional systems such as U2 flights. To 
make such comparisons we must be able to evaluate the effectiveness of 
each type of surveillence system. This paper will be concerned with eval- 
uating satellite surveillence systems. 

What elements are important in the evaluation of a satellite surveil- 
lence system? This is a difficult question to answer. The appropriate 
elements must include both the nature of the phenomena to be observed, and 
the available hardware, including the satellite, its components, and ground 
support equipment. The nature of the phenomena to be observed may be such 
that one satellite of a particular type is sufficient. Or, it may be 
found that multiple satellites of one or more types are required. The 
appropriate number will depend upon the effectiveness of each individual 
satellite. Thus, an initial appraisal of a system should be based on an 
individual satellite's performance. 

What features of an individual satellite's performance should be ex- 
amined? Two features appear to be dominant. For those phenomena which 
occur randomly in time and location on the surface of the earth, it seems 
desirable to make some statement about the expected proportion of the total 
phenomena of interest the satellite will observe. For other phenomena we 
are more interested in seeing the result of an event rather than the event 
itself. Thus, we would like to know the average elapsed time before a 
satellite observes a location which has suddenly become of interest. (For 
example, if on the average 100 natural disasters occur in a year in the 
area covered by the satellite, then we may want to know the average time 


lapse between occurrence of a disaster and observance of the place of 


disaster.) 

The present procedure for satellite system evaluation is to construct 
an elaborate model describing the location of the satellite at all times 
for a given orbital inclination, altitude, etc. This approach requires 
that every combination of orbital inclinations, sweep width, altitude, 


etc., must be investigated to have a complete data bank for estimating 





the frequency of observation or the time before an event is observed. 
This, by nature of the model's complexity, requires the use of a high 
speed computer and much computation time. For example, if each of three 
parameters is to be investigated for 10 different values, then 1000 separ- 
ate replications would be required. 

The author believes that an alternative approach, by using a simpli- 
fied analytical model, may prove to be just as useful for most appraisals 
of satellite surveillence systems. Appropriate analytical expressions, 
even if approximations to the actual situation, usually help to point out 
the important parameters, variables, etc. AS a consequence, a Simplified 
model like that discussed in this paper may point out some possible con- 
clusions not immediately apparent in a computerized model. For example, 
it was mentioned that the United States launched seven surveillence satel- 
lites during 1965 which were characterized by polar orbits. It will be 
obvious from the model to be presented in this paper that for reconnaissance 
purposes a polar orbit is extremely inefficient unless the polar regions 
themselves are of particular interest. 

Basic to the approach is the development of a fairly simple analy- 
tical expression for the likelihood of seeing any event of interest. It 
would be convenient if we could find the probability that the satellite 


passes over a chosen ground position during any randomly selected orbit-- 


but since the satellite's position is completely deterministic we can not 
speak of the "probability" that the satellite passes over such a position. 
(For a selected orbit either the satellite passes over the position or it 
doesn't.) However it does make sense to speak of the "fraction of orbits" 
in which the Satellite passes over a selected ground position. This dis- 
tinction between "fraction" and "probability" is a fine one, but an impor- 
tant one. For example, if a satellite sees London every seventh orbit then 
it sees London on one-seventh of the orbits, but the probability it sees 
London on any randomly selected orbit is not one-seventh. 

The likelihood function to be derived will, then, describe the frac- 
tion of orbits that the satellite spends in "positive" orbits, where a 
positive orbit is defined to be one such that the satellite passes over 
the ground position in question on that orbit. If, for example, the ground 
position is Beirut, then the fraction of positive orbits would be the ratio 
of the number of orbits in which Beirut is within the sensors sweep width 
to the total number of orbits (per some time interval, e.g. the life of 
the satellite). If P represents this fraction, then, using reasonable 
assumptions, P can be shown to be a function of the satellite's orbital 
inclination, the sweep width of the sensor, and the latitude of the point 
of interest. P is also a function of the sSatellite's altitude, in that 


the sweep width is usually a function of the satellite altitude. 


oe Formulation of the Likelihood Model 


The formulation of the likelihood model is facilitated by the 





following sketch. 











Satellite motion in general is confined to being great circle mo- 
tion=. [hat is, the ground track of the satellite describes a great 
circle on a stationary globe. Controllable parameters are altitude, 
orbital inclination, and sweep width. 

ax surveillence purposes, upper and lower bounds exist on satel- 


lite altitude. The lower bound is caused by the earth's atmosphere. A 


Satellite must be outside the atmosphere to remain in orbit. The upper 

limit is dependent on the resolution capabilities of the satellite sen- 

sors. The lower and upper bounds are presently about 100 and 300 miles 

respectively. The 300 mile maximum altitude means that the time for one 
orbit will be less than 1.4 hours. That is, while the earth rotates one 
revolution the satellite has orbitted about 17 times. This formulation 

will assume that altitude and sweep width are directly related and that 

knowing Sweep width is equivalent to knowing satellite altitude. There- 
fore no further mention of altitude need be made. 

The orbital inclination angle © is equal to the extreme latitude 
attained by the satellite. Once launched, we assume @ remains unchanged 
for the life of the satellite. 

The sweep width 2r is the width of the path on the earth which the 
Satellite sensor Sweeps out as it proceeds on its orbit. For photographic 
sensors 2r would be the width of the area photographed on one exposure, 
for radar it would be the length of one scan, scans being perpendicular 
to the direction of the satellite motion. 

We let & represent the latitude at which an event occurs, and since 
we have no control over the location (or time) of events of interest, « 
is an uncontrolled parameter. The ground position at latitude & which 
is of interest will be designated as A, where A is either selected ran- 
domly or non-randomly, depending on the context. 

General Assumptions of the Model 

fhe following general assumptions are basic to the likelihood model: 

(i) The satellite under consideration remains in orbit for long 

periods of time, upwards of 1000 orbits. The life of a sat- 


ellite will be denoted by L. 
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(ii) An event of interest is equally likely to occur at A at any 
time. Given an event occurs during the time L, then we 
assume the random variable denoting the time of occurrence 
is uniform on (0,L). 

(iii) 9 cannot be varied after launch. 

(iv) The orbit is circular. 

(v) Every point on earth which lies within the sweep width of the 
satellite's sensor as the satellite passes overhead is seen 
with probability 1. 

(vi) The earth is a perfect sphere. 

(vii) The sensor points normal to the earth's surface at all times, 
and can see a distance r to each side of the satellite's 
ground track. This varies with the actual sensor aspect in 
designed systems (in general the sensors are aimed to the 
side) but simplifies the trigonometry, without changing any 
basic concepts of the systen. 

(viii) Sensors are either photographic or radar. In either case, 
the ensuing model assumes perfect visibility. (Brief men- 
tion is made in Section 4 of how one might incorporate the 
weather problem, and consequently the problem of general 
atmospheric attenuation, into the model) 

The Question of Accuracy 
The formulation to be presented does not yield extreme accuracy of 
results. Extreme accuracy seems unnecessary in the problem of initial 
system evaluation. For example, if we determine, through the model, that 
the mean time to next sighting of a ground position A is 2 hours when in 


actual fact it is 2.078 hours, then this probably would be of small con- 
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sequence in any decision made regarding the system. Certainly it would 
be misleading to present results implying accurately to six or seven fig- 
ures when basic assumptions preclude accuracy to more than two or three 
figures. Computer formulations specifically designed for the surveillence 
evaluation problem may claim accuracy to six figures. [3] These figures 
may or may not be that accurate when compared to an actual satellite in 
orbit, but any initial evaluation of a satellite system would probably have 
remained the same had the figures been rounded to two significant figures, 
which is approximately the accuracy inherent in the proposed model. 
The Likelihood Function 

We will denote the "fraction of positive orbits" by P(@,a« , 1). 
However P(@, om, r) will often be referred to as the likelihood function 
P. As the notation implies, P depends on three parameters: orbital in- 
clination 6, latitude @& of the ground position A to be observed, and the 
sweep width 2r. 

The analytical expressions for P(@,« , r) are derived in Appendix 
I and are given below. The derivation is based on the assumption that P 
is equivalent to the ratio of the sweep width measured along a latitude 
to one-half the circumference of the latitude. (i.e. TS in the previous 


LA 
sketch) 
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where § = =: and R is the radius of the earth. 


Figure 1 illustrates the behavior of P as a function of & for 
several values of 6 when r = 150 nautical miles. The figure shows 
that P attains its maximum when & = @-8. This fact can be con- 


firmed by an investigation of (1). 
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3. Applications. 
Mean of Time Until Next Sighting 

One of the present parameters used for appraising a satellite sur- 
veillence system is the mean time between sightings. + From our model the 
mean time between sightings is simply — . However if location A we 
suddenly become of interest we probably would not care about the time 
lapse between the last sighting of A and the next one, but rather on the 
time lapse from the present until the next sighting. These two parameters 
are by no means equivalent. 

In this section we will develope several estimates for bec mean time 
from a given instant until the next sighting of a specified point A. In 
these developments, we will assume that phenomena of interest are equally 
likely to occur at A at any time during the life of the satellite. 

If Ty is the time lag between occurrence of an event at location A 
and the first sighting of A after the occurrence, then we are looking for 
the mean value of T,, denoted by M(T;,)-. An exact expression for M(T,) is 
difficult to obtain, but various estimates can be made. 


To see why the exact expression is difficult to obtain analytically 


consider the following sketch. 


lror an example, see [1]. 
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If we consider an upward (downward) pass through the latitude @ to be one 
such that the satellite is moving from & to a more northerly (southerl]ly) 
latitude, then the satellite can see point A on any given orbit in one of 
two ways. Either it sees A on an waged pass or on a downward pass. The 
time between an upward and downward pass sergbusdle on the latitude &. If 
T' deontes this time and T is the time for a satellite to make one con- 
plete orbit, then the discrete points in time at which sightings can pos- 


sibly be made occur in a sequence Ty» Toseeeces te » where the times 


between Sightings are ‘Tae To 


-_ ie ‘ies Ty se To, eee 
In reality, sightings do not occur in a neat sequence Ty» Too-++s Ty» 
-..- aS described above. That sequence merely tells when sightings are pos- 


sible. The actual sightings may occur in sequences which are exceedingly 





difficult to describe analytically. For example, as indicated by the 
previous sketch, a sighting occurs when on an upward pass the satellite 
passes over point P while A lies between Q and Q'. Since the distance 

QQ' is on the order of several hundred miles, the speed of the earth's 
rotation makes it possible for A to transit QQ' in much less time than it 
takes the satellite to make a complete orbit. Thus the satellite may have 
just crossed & while A was just west of Q, and A will pass Q' before the 
Satellite makes another upward pass. A "sighting possibility" would thus 
be missed. 

Rather than attempting the difficult task of predicting the actual 
sequence of positive orbits, we will develope estimates for M(T,). If 
tests show these estimates to be reasonably accurate, then they should be 
almost as useful to the decision maker as the true value. Let us consider 
three possible estimates. 


Estimate I: If we assume actual sightings occur at equally spaced inter- 


16 


vals of time, then the spacing can be obtained from the likelihood function 
P. The assumption that events of interest occur with equal likelihood at 


any time during the life of the satellite implies that if an event occurs 








during the time CT; Tj4 )) then it is equally likely to occur at any in- 


stant in that interval; that is, the conditional time of occurrence T is 





uniform on (T;; T ).- Since the length of this interval is ; MCT, ) 5 


i+] 


an estimate for M(T,) > is 


Intuition suggests that this may be a good estimate near the equator 
where the time between an upward and downward pass, and.therefore the time 
between possible sightings, is the constant “1. Since the times of pos- 
sible sightings are equally spaced, the time between actual sightings may 
tend to be equally spaced also. But it seems equally intiutive that as A 
moves away from the equator the estimate will become less reasonable. 
Estimate 2: Arbitrarily define a probability model for the orbiting sat- 
ellite, wherein P, instead of being the fraction of positive orbits, is 
defined to be the probability of a positive orbit. Further assume that 
all orbits are independent of each other. This model is then equivalent 
to one with Bernoulli trials with probability of success equal to P. Fur- 
thermore the number of trials to first success for such Bernoulli trials 
is a geometric random variable with mean —— where P is the probability of 
success on each trial. Equivalently, in terms of our model, the number of 


orbits to the first success (positive orbit) is _! and therefore the 


estimate for the time to first positive orbit would be 2 





wen 


Intuitién suggests that this estimate should be good when «& ap- 
proaches 6, but not so good near the equator. The time between an upward 


and downward paSs approaches zero as &->©@, and between downward and up- 
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ward approaches T,- Thus the time till possible sighting will be twice as 
long as that at the equator. 

Estimate 3: Since neither of the above estimates seems reasonable for all 
values of am , we return to the physics involved to see if an assumption 
can be made which provides awbetter estimate. 

If a ground position A is at latitude q@ we have said it can be seen 
during that orbit at one of two times, it can be seen during the upward 
pass or the downward pass. It is relatively easy to obtain an exact ex- 
pression for T', the time between the upward and downward passes. Knowing 
T, and T we can describe the sequence of times at which sightings are 
possible. We also Know that the fraction of positive orbits is P. 

Using these two pieces of information we assume that actual sightings 
occur in a sequence such that the time between sightings is proportional 
to the time between possible sightings, where the constant of proportion- 


ality is such that the fraction of positive orbits is P. We get (see 


Appendix IT) 
M (7%) = ao {leos"/(cof G fan «)]* =f 


[77 - <os(cot @ tana)]} (3) 





A T 
Notice that for & = 0,M(T,) becomes = which is identical to 





estimate 1, and for & = 6, M(T;) becomes ~ which is identical to 
estimate 2. Thus estimate 3 varies with e& and also fits what inttition 
suggests should be the case at the extremes of the range of & . Further 
evidence of this estimate's reasonability is given by Appendix III where 
MCT) is compared with results obtained from a computerized simulation 


developed by the Naval Missile Center. (3] 
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Location and Number of Readout 
Stations Required for a Specified Level of Contact 

Any study of the effectiveness of a satellite surveillence system 
must include an appraisal of problems of transferring ‘the information 
from satellite to evaluation centers located on earth. To consider a sat- 
ellite which orbits for long periods of time we should discuss readout 
stations. Readout stations are radio receiving stations capable of re- 
ceiving signals from a satellite whenever the satellite is within a cer- 
tain distance of the readout station. We will refer to the circular area 
wherein such reception is possible as the transmission area. As the sat- 
ellite passes through a transmissidn area it transmits information accumu- 
lated since the last pass through a transmission area. The author knows 
of no present study which analyzes the optimum location of such stations 
or the number needed to ensure adequate transmission time, therefore it 
seems worthwhile to consider this question. 

P(@,a , r) can be used to determine the percentage of time a satel- 
lite will be in contact with a particular readout station. The procedure 
we will use is: 

1. Find the diameter of the transmission circle within which 
radio signals from the satellite are receivable. 

2. Estimate M(t), the mean time a satellite takes to cross the 
transmission circle given it intersects the circle. 

3. Determine the value P(@, & , r), where & is the latitude 
of the transmission circle and r is the radius of the transmission circle. 

4. Estimate the mean time per orbit available for transmitting, 
by P\ M(t). 


Knowing the rate of transmission of information, a reasonable esti- 


ng 


mate for the number of readout stations needed to maintain a specified 
level of information transmission can be obtained. 
To find the diameter of the transmission circle consider the fol- 


lowing sketch. a 


In this sketch h is the height of the satellite above the earth's sur- 
face, and a is the minimum reception angle above the horizon. 


From the law of Sines, 


Sin (angle ABC) a sin (ar+Z) =. COS a 





R R Ya V7 Ar } 
T here tore , Angle ABC = sin (- 4 = =) 
and , | i “ oes ae 


Thus, the diameter of the transmission circle, d, is 


d= Z2(R+H) eo , (4) 


—_— 


Thus we see that a readout station can receive signals from the satel- 
lite if the satellite is within a circle with diameter 2(R + hp with 
center at the readout station. If we not let — the expression for 
P(O, %, r) equal 1(R + hp » then we will have the expression for the 
fraction of orbits wherein the satellite passes within the transmission 
circle. 

To find the mean time available for transmitting given the satellite 
passes through a transmission area, we will assume that the distribution 
of paths which the satellite makes as it crosses the transmission area is 
uniform; that is, the satellite is equally likely to cross the circle on 
any chord. Finding the mean length of path based on this assumption is not 
difficult. If the radius of the transmission circle is r then the mean 


distance travelled while in the circle is approximately . (See Appen- 





wer 
2 
dix IV). 

From the above results, we can obtain an estimate for the mean time 
per orbit spent in the transmission area of a single readout station. It 
is 


M(t) = -- P(®, M%, rr); r= (R+H) 





V 


ground velocity of satellite. 

Because P is maximized when X= @-§ and M(t) is linear in P, we 
can maximize the mean transmission time per orbit by placing the readout 
station at latitude @-§ . 

We next consider the question of how many readout stations &re re- 
quired to attain at least some specified mean transmission time per orbit. 
We can reaSonably assume that all stations would be placed close to lat- 
itude @-6§& . we also assume the stations are spread out on that lati- 


tude to avoid a satellite passing through two stations in rapid succession 


21 


and then not passing through another for many orbits. Under these assump- 
tions a reasonable estimate for the number of readout stations N required 


to average t’ minutes of transmission per orbit is 


t 


N= —< 
M(t) 


If a value of t’ is designated as the minimum mean time per orbit needed 
in order to receive some specified fraction of the total information ac- 
quired by the satellite, then N tells us how many readout stations are 
required. Such information would be needed to estimate the cost of a pro- 
posed satellite system. 
Estimating the Number of Intersections of a Transitting Ship 
by a Satellite 
To find an estimate for the number of intersections of a moving point 

on earth, consider a ship going from point a to destination b, where a 
and b are at latitudes A and 8 respectively. Let us assume: 

1. The ship proceeds at constant speed s. 

2. The ship moves on great circle path. 

3. The distance d from a to b on this path is known. 

4. The track of the ship does not leave the area of satellite 
coverage. 
First, to avoid complicated expressions derived from differential equa- 
tions giving the instantaneous latitude of the ship as it proceeds, we 
elect to use a constant, 0 , to represent the time averaged mean latitude 
of the ship; that is, 0 is an estimate for the latitude which, if used to 
replace & in P(®,& , r), would give the same result as if we continu- 
ously evaluated P at each instant of time during the ship's transit. 6 is 


defined as 


Ze 





Since the ship's track is a great circle track and the ship's speed is 
constant, this is a lower bound for the ship's time averaged mean latitude. 
A conservative estimate for the number of intersections will consequently 
be obtained. 

The derivation of the estimate of the total number of sightings for 


a ship transitting from a to b under the above assumptions is the follow- 


ing: 
Let h = satellite altitude 
t = ship's transit time from a to b 
v = satellite's speed 
g = gravitational constant. 
Then a 
v= =e »5 and t = ae 
new ao 20m (R + h) 
Oo Vv 


and thus N, the number of orbits during the transit is 


N = ° 





If I is the number of intersections during the ship's transit then an 
estimate for I is 

A 

I = NP(@, DS , r)% 
Example: If the ship goes 1000 miles from 10° to 20° at 12 knots, the 
average number of sightings by a 200 mile high satellite with orbit 
@ = 50° and r = 150 will be 


IT=P-N=NP(O,D , vr) = (63.2)(.035) = 2.2. 
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4, Areas for Further Study. 
Having derived an expression for the fraction of positive orbits one 
can investigate various areas other than those already mentioned. The 
following are areas which seem to be of obvious interest in the evaluation 
task. 
1. Finding the optimum orbit based on some measure of effectiveness 
is a problem of considerable interest. Figure 2 shows graphs indicating 
the distribution of major nautical activities as a function of latitude. 
Further work could be accomplished correlating P and these distributions 
to optimize the orbit. For example, if the measure of effectiveness is 
‘number of contacts per unit time'' then we could proceed as follows. Des- 
cribe a graph of figure 5 by a function f(« ), where f is either a discrete 
or continuous function approximating the graph. If we consider southern 
latitudes to be negative, then we can let 
a | o 
Fla) = {fda | or F(a) = 3) Ha), 
= 


77 


7 2 
gz 


depending on which is applicable. Then, let 
G(e ,r) = (PCG, 4,7) d Fla), (10) 
of 


If G is the maximum of G with respect to 9,, then that 8; which yields 
G is the optimal orbital inclination for our chosen measure of effective- 
ness. 


The distribution of naval activities with longitude may be combined 
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Figure-2. Distribution of Various Naval and Maritime Activities 
with Respect to Latitude. [2} 
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with the distribution for latitude. The optimum orbit for observing speci- 
fied regions may then be quite different than if latitude distribution 
only are considered. 

2. In studying the effectiveness of a satellite as a weapon delivery 
vehicle, we need only replace the half sweep width r by the lethal range 
r' of the weapon carried by the satellite. Then, if the ground position 
of the satellite is within the lethal range r' of the target, the weapon 
could be dropped. P and M(T, ) would have the same form as above, but 
M(T],) would be the mean time until a weapon could be dropped on a target 
after the decision to destroy the target had been made. One might cor- 
relate this with the possibility of inflight satellite course changes. 

3. Investigating countermeasure techniques may prove of interest. In 
this respect figure 1 gives some indication of P and its sensitivity to 
latitude in the extreme regions near 8. Knowing the inclination of an 
enemy's orbiting satellite and correlating that information with the impli- 
cation of figure 1 would indicate ship routing procedures we could use to 
avoid detection. We obviously would try to route our ships to Stay clear 
of the latitudes near @-§ , where P is largest. 

4. When data is available on global cloud cover distributions, P can 
be modified by the probability of cloud cover and/or darkness obscuring 
the sensor's viSion (as in the case of photographic sensors). As a con- 
sequence, P would become closer to being a probability than a fraction. 

As a result of cloud cover being a function of longitude as well as lat- 
itude, P would also have to be dependent on longitude. The problem of 
predicting cloud cover is obviously complex and any model developed to 
describe the distribution of cloud cover would probably have to be a con- 


puter simulation. 
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APPENDIX I 
DERIVATION OF THE FRACTION OF POSITIVE ORBITS 

The function P(@&, & , r) is defined to be the fraction of 
positive orbits to total orbits given a satellite in orbital inclina- 
tion 8, a sensor sweep width 2r, and a ground position ae Sa . 
The basic rationale for finding P is to consider the circle on the 
sphere made by latitude & , and to find the ratio of the path width of 
the satellite measured along latitude & to one-half the total acerca 
Then if the ground position A is on latitude & , the fraction of times 
A is seen will be proportional to this ratio. 

To obtain an analytic expression for P based on this rationale 
‘ consider the following sketch, in which capital letters denote points 
on the sphere, lower case letters wenste angles of spherical triangles, 


and Greek letters denote sides of spherical triangles, a side being the 


usual radian measure of the central angle subtended by that side. 
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If we could determine the length of BB', then divide BB' by FG, (one 
half the circumference of the circle made by latitude @ ) we would have P. 
However it is not necessary to find BB'. We get the required results by 
finding, instead, angles a and c, subtracting these from tr to get angle 
b, then dividing b by TT. 
To find angle a we will make use of the fact that triangles AND and 
CBD are spherical right triangles. Furthermore, side Y of AND is 1 - 6. 
Thus if we can find side A then Napier's Rules for spherical right tri- 
angles will give us angle a. A similar approach is used to find angle c. 
We first review Napier's Rules. Given a spherical right triangle as 
shown in the sketch below, where 4, > and Y represent the central 
angles subtended by the respective sides of the triangle, then the rules 
can be summarized as follows: 
Take the five parts, excluding the right angle, and con- 
sider them to be in the circular arrangement shown. If 
"co-" means "complement of" then: 
(1) the sine of the middle part equals the product of 
the tangents of the adjacent parts; 


(2) the sine of the middle part equals the product of 


the cosine of the opposite parts. 
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To find angle a we first note that g is a right angle since A is 
the northernmost point of the orbit, and h is a right angle by con- 
struction. The sweep width is 2r and the central angle subtended by r 
ie then ‘ ea 
From Napier's Rules, 


sin § = cos (co-§) cos (co-e) 


and therefore 


sin § 
sing = S22. (11) 
Similarly, 
sit7 P 
7 (12) 
But ; 
+5 = = ae . . (/3) 
2 
From (£1), ({2), and (13), m 
sinP sing _ sin¥ cos (x +) (14) 
sin Xr = a ag _ s/77 & : 


After some manipulation of (14), the expression for A is 


A “ff sin § + cos (§4+A) sin 
A= cot Me Tay So ee ° (15) 





Again using Napier's Rule, 


cos 2 = tan cof A . (/6) 


From ({5) and ({6), 


ie sin §& + cos (§tA) sin 
Qs Cos ‘| sec pf eee ee PI]; (17) 
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To find angle c consider the following sketch. 


N 





) . y 


Due to symmetry between northern and southern hemispheres, B'G equals 


B''G''; therefore angle c equals angle k. But angle k obviously equals 
angle f; therefore to find c we need only find f. Now f can be found in 
the same way we found angle a by merely using spherical right triangles 
AND'' and C''B''D*' instead of triangles AND and CBD. Therefore we aan 
use feiien (17) directly by replacing (A +§ ) with ( X»' + é*) a Bue 
(r' + §') must equal T-( A +§ ) because (€* + XX") equals (7 -Y ) 
and > equals ( A + §€ ) by symmetry. Therefore, after the appropriate 


changes in (17), we get 


oF sin & & Cos (1r-§-A) sin? sj: (18) 
aa °°" [see ¥[- sin (7 -§-4) -/ 


We are now ready to obtain an expression for P(@,a&,r). We 


said P would be equal to — » where b 1s obtained from 
b= w-(a+ec). (1%) 


Using the relationships 


E+A 


if 
Ny 
j 
R 


; ie 
3 a7 9) 


and PCOA%X, 41) = £ 


in equations (17), (18), and (19), we obtain / 


} =) gse2 (57 + SINK COSC 
PCO, 4,1) = - pfs COS o si & me Se ») + 





On a reexamination of the derivation we realize that (20) is valid only 
for M&S @-$ . Because angle a is zero for xe> & -§ , (20) must 
be modified by dropping the first term in ¢-? for O-$S%S Or §, 


Finally. P must be zero for am>&+§ . Therefore the expression for 








P over the range 
(0 S + SINK COS @)|* 
O< OsF Me 


/- 3 5 [orl 4 
~t “Ss fp; ?94u<O-5 
— 


Poe. 


Ce 


The behavior of P over the range 0 £<& “eZ is illustrated in 


figure i. 


APPENDIX II 


Derivation of Estimate 3 for Mean Time to Next Sight ing 


Consider the following sketch: 


N 


We will consider our line of sight to be always across the plane 


of satellite rotation, and let the earth turn relative to that plane. 


let D 


D' 


Ty 


= time for’single revolution of earth in hours. 

= time in ee earth to turn far enough so that A (fixed 
on earth) moves twice the distance QR. (Thus D' is the 
time between a possible Sighting on an upward pass, and a 
Beeeivie Sighting on a downward pass.) 


= ‘Orbital period in hours 


From the sketch, 


D 


= arc QR_ p = _angle g p 
arc PR wT 


But, by Napier's Rule, 
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cos g = cot @ tana. 


-1 
Therefore D'= scos7= (cot_Y tana ) D. 


Tr 

We want to find an estimate for the mean time to next sighting of 
location A, given A has suddenly become of interest... We are unable, with- 
out great difficulty, to describe the sequence of actual sighting times. 
Therefore we would like to define a sequence of times which we can use to 
approximate the sequence of actual sightings, and upon which we can base a 
reasonable estimate for M(T}). 

We will assume that the sequence selected should satisfy two condi- 
tions. 

Condition 1. The sequence must satisfy the condition that the 


mean number of orbits between positive orbits inn, where n = : 


Condition 2. The sequence should be proportional to the sequence 
of possible sighting times. 
It is obvious why condition 1 is desired. Condition 2 is specified be- 
cause it will require the sequence of assumed actual sightings to be 
spaced in a way which suits our intuition. Thus, if the possible sight- 
ings occur in a sequence such that the time between possible sightings 
has the sequence 

Menee 25 ir, | hr, 23 Rts. + «eae, 
then we would intuitively expect actual sightings to be spaced such that 
usualip tao sight ings would occur relatively close together, followed by 
a long period with no sighting, then two more sightings relatively close 
together, etc. 

Therefore if the times between possible sightings has the sequence 


B= D', D-D', Ds D-D' ,...- 
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then we want a constant K such that the sequence KZ satisfies our two con- 
ditions. There is obviously only one sequence which satisfies these two 


l Therefore our next task is to find that se- 


conditions simultaneously. 
quence. 

To find the sequence KZ we need only find a K which, when multiply- 
ing Z, satisfies condition 1, because any Sequence KZ where K is constant 
will automatically satisfy condition 2. Now the sequence Z "cycles" on 
every second sighting possibility with a cycle time of D and therefore 
there are two possible sightings in a cycle of Z. We would like kZ to 
"cycle" on every second sighting also, but the cycle time must be equal to 
2nT5 >» Since there are to be 2n orbits per every two sightings. Therefore 


2nT 
to get KZ we must multiply Z by = Thus: 


i= al 
D 





We then have 





i 
AZ= f £22 p Ml po) a a ae 
Eo, Bhiono), BE p; Bah}. 


To show that this sequence satisfies condition 1 we consider a single 
cycle (22k p- Lk (p-D’) 2 The mean time between sightings is, 
e D 


then, 





/ Lat 70 Ay 
s[— b+ 2% (p-0) =, 


and therefore the mean number of orbits between positive orbits is n. 


lie have fixed the mean time between points in the sequence and 
the "cycle" of the sequence. 


36 





Having defined KZ, the sequence of times between assumed actual 


. ° te ° : * 
Sightings, we can say that assumed positive orbits occur at times Ty» 


x * 
To) Tarceeece where 


7k To) + Rn 2) 3 € odd 
cf = 

* D-2D : C ver) 
Ey + 2nb( 22 a 


he. f7* hl ae TD oe 





where T. indicates the time of the jth assumed sighting and the time be- 


C4. 


tween T; and noe is defined to be one cycle. Then assuming the time of 
occurrence of an event of interest at A is uniform over the satellite's 
life, we can-say that the event of interest is equally likely to occur 
; ' , kk 

at any time T contained in (T;, Ti42)° 

We are trying to estimate the mean time to next sighting given an 
event of interest has occurred at A. However, since each cycle is iden-. 
tical, we can equivalently find the mean time to next sighting given 


rT é (Tr, igo) The density of T, where T is the random variable denoting 


ay 


the time of occurrence of the event of interest given T € ies Teas 





is then 
/ 
t)= ——. 
f(t) Zn 
Therefore 
M(t) = Te, time trom T to Te, 
ee ae a Rew =; 


7g, [mean time trom T to i 
Pha weiter 7 é he; 7 


Le/ 2 


7, |y* * x), B2_[-* 
= ee [i I,h-E(r/Te (et ish) se Tena -E(TITE(ES, Tea) 





- 2 ITS - Fe Pi " 7&2 [2 - fea Tit | 
Py: Lt/ 2 2rniz LTA ae ° 


Substitution of the appropriate expressions for the T3"s in (22)- results 


2 Lo)» (o-0)'] 


DP 


in 





M (i) = 


and, upon substitution for D’' from (21), 


A O b, « 
MT) = 22 [2s ‘(cot 6 tana) oe 


fer — cos! (cot P tan o7*/ (25) 


Thus, we obtain a general estimate for M(T;) as a function of @,am , and 


Te 
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It is of interest to note that when & = 0, (23) reduces to 


A ube 
M(T, ) = pee 


and when & =@, (23) reduces to 


A fis 
M(T,) = 





These limiting forms are identical to estimates 1 and 2 respectively, 
which is a desirable result, as estimates 1 and 2 seem appealing at these 


limits. 
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APPENDIX III 


Comparison of Theoretical and Numerical Estimates for M(T}{) 


Our analytical results are compared with simulated results for the 


following data: 9 = 50°; r = 150n.ms; B® = 1.6@ms; 05 orbure, 


fe) 
ox Analytical Results Simulated Results 
0 32 hrs 35 hrs 
20 2.2. eles 22% Tirs 
40 16.5 hrs 16.7 hrs 
50 12.5@hrs 10.5 hrs 


The simulated results were obtained from a Stanford Research Insti- 
tute study which computed the times of intersection for each of some 1800 
fixed locations. [2] To obtain M(T,) from these times of intersection 
for each location the occurrence of an event was assumed to be uniformly 
distributed over time during the 105 orbits considered. M(T,) for loca- 


tion j was then found by the formula 





where XY) lumalia, °°) are the times of intersection, Aj is the total 


n, 
time considered for location j(F3 = time between the first and last in- 
tersections for j). This formula is easily verified as a formula for the 
time to sighting of an event which occurs uniformly in (0, 5) when times 
of sighting occur at times W,,&,) --- Xn; » Then the overall value of 

MCT, ) was found by averaging the M(T). Since locations were scattered 
throughout various latitudes only those within a small band were considered 


to find M(T,) for that band. For example, all locations between 39° “anid 


41° were used to find M(T, ) for 40° latitude. In the case where two cities 


40 





were so close together that considering both would create bias, (Los 
Angeles and Long Beach) only one was used. 

The results are extremely close with the notable exception of the 
range where we would expect our estimate to yield the best results, i.e. 
at & = @ = 50°. ‘The difference here is most likely due to the fact 
that’ the model in this paper assumes the satellite sensor "looks straight 
down" whereas the computer model has the sensor's field of view skewed 
toward the north side of the path. This difference between our model 


and the computer model can be seen from the following sketch. 





skewed sensor 


S 7 minetrie Selrsor 


The satellite in the latter case, when it is over latitude 50°, 
could have its sensor centered on a latitude as high as 54°. While a 
four degree error is negligible at moderate latitudes, as & approaches 
6 the difference becomes quite large. If one fraction of positive orbits 
from the computer simulation is computed and used in our formula instead 
of PCO, x; r), then M(T,) for * = 50° is found to be 10.7 hrs. There- 
fore the variation between our results and those of the computer simulation 
% appears. to be caused by the difference in fraction of positive orbits due 


to slightly different physical modeis. 


APPENDIX IV 


Derivation of Mean DiStance to Cross Transmission Area 


We want to find the mean distance X that a satellite travels when 
it transits a transmission area of a ground readout Station. We can 
estimate this distance by assuming the transmission area is a plane cir- 
cle, and that the track of the satellite is a chord through the circle 
.with the distance from the chord to the center being uniformly distribu- 


ted on (0,r). To facilitate the discussion, consider the following 


sketch. 





In terms of the sketch we want to find E(X), the mean of X, when 


the random variable Y is uniform on AC, that is, 
; of yer 


l 
cao) je = 
f(y O ; otherwise . 
From the sketch we get 


xX’ = (re-v)* 
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and, because Y is uniformly distributed, the random variable X' has 


the distribution function: 


Eade 1- BY * (PR!) 


2 ry2\Z 
ja (208) 5 Of X'EF. (24) 


But X = 2 therefore the distribution function of the random variable 


X is 


_ (X)*)\4 | 
FE, (#) = /- (cB) P ‘ O 2 Kiaeers (25) 


- It is then easily shown that 


E00 = EF " 


Equation (26) was obtained assuming a plane circle, whereas the actual 
transmission area is a segment of a sphere. The result is a reasonable 


approximation, however, because the transmission areas are small enough 


that the curvature of the earth is negligible. 
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